Quasiparticle bandstructure effects on the d hole hfetimes of copper within the GW 

approximation 
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We investigate the lifetime of d holes in copper within a first-principle GW approximation. At 
the Go Wo level the lifetime of the topmost d bands are in agreement with the experimental results 
and are four times smaller than those obtained in the "on-shell" calculations commonly used in 
literature. The theoretical lifetimes and bandstructure, however, worsen when further iterative 
steps of self-consistency are included in the calculation, pointing to a delicate interplay between 
self-consistency and the inclusion of vertex corrections. We show that the Go Wo success in the 
lifetimes calculation is due to the opening of new "intraband" decay channels that disappear at 
self-consistency. 

PACS numbers: 71.15.-m, 71.20.-b, 79.60.-i 



Very recent experimental and theoretical results on the 
quasiparticle lifetimes in noble and simple met- 

als show that our present understanding of the electron 
dynamics in real solids is far from being complete [Q. 
For electrons above the Fermi level (hot electrons) time- 
resolved two-photon photoemission experiments show a 
non quadratic behavior of the lifetime, in contrast with 
the Fermi liquid theory prediction [||. For occupied 
states, the number of possible scattering events increases 
rapidly as their energies decrease below the Fermi level, 
and the agreement between photoemission experiments 
and theoretical results worsens |^]. An open question is 
whether this discrepancy comes from effects beyond the 
approximation used in the calculations (i.e., beyond GW, 
where vertex corrections are neglected) , or it comes from 
the way used to solve the quasiparticle equation for a 
given self-energy. In the present work we address quan- 
titatively the latter point, the inclusion of vertex correc- 
tions being beyond the scope of this paper. 

Density Functional Theory [Q (DFT) has become the 
state-of-the-art approach to study ground state proper- 
ties of a large number of systems, going from molecules, 
to surfaces, to complex solids The success of DFT 
is based on the idea that the spatial density of a sys- 
tem of interacting particles can be exactly described by 
a non-interacting gas of Kohn-Sham (KS) independent 
particles, moving under the action of an effective po- 
tential which includes the exchange-correlation potential 
Vxc- Compared with experiment, the usual local den- 
sity (LDA) 1^ approximation to the DFT yields semi- 
conductor bandstructure which systematically underesti- 
mate the band gap, while for noble metals the discrepan- 
cies are both k-point and band dependent jlOj. Another 
important drawback of the use of DFT eigenvalues as ex- 
citations (bandstructure) energies is that they are by con- 
struction real; no lifetime effects are included. An alter- 



native approach is Time Dependent DFT where all neu- 
tral excitations are, in principle, exactly described [pT| . 
However, quasiparticle lifetimes have not been considered 
so far within this approach. 

Many-body perturbation theory allows one to obtain 
band energies and lifetimes in a rigorous way, i.e. as the 
poles of the one-particle Green's function G . Those 
are determined by the solution of a Dyson-like equation 
of the form ill: 
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containing the non-local, generally complex and non- 
hermitian, frequency dependent self-energy operator E. 
The poles of G are the quasiparticle (QP) energies e^jf , 
that from Eq. (|l|) correspond to the generally complex so- 
lutions of the equation: = £'„k (^nif) • "^^^ ""^^^ P^'^^ 

(DP. 

of e^j^ gives the quasiparticle bandstructure; whereas 
the imaginary part yields the inverse quasiparticle life- 
time. In the present work E is evaluated according to 
the so-called GW approximation, derived by Hedin in 
1965 jl^ , which is based on an expansion in terms of 
the dynamically screened Coulomb interaction W (cu). In 
the first iteration Go and Wq from DFT calculations are 
used to compute E as Eq = iGgWo. Unlike semicon- 
ductors, the case of noble metals has been studied only 
recently ||, |l^, |l^ . The lifetime of hot-electrons in cop- 
per has been calculated with the GW self-energy eval- 
uated at the DFT zero-order energies ^ (namely "on 
mass-shell" approximation). This approach is based on 
the assumption of vanishing QP corrections of the DFT 
states while, very recently, large QP effects on the oc- 
cupied bands of copper have been found 0. The on 
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mass-shell" approach applied to the hole lifetimes in- 
dicates that d-holes in copper exhibit a longer lifetime 
than excited s/p electrons. The quantitative comparison 
with experiment however, has shown a large overesti- 
mation of the experimental lifetimes measured by means 
of photoemission spectroscopy. In this paper we go be- 
yond the "on mass-shell" approximation, and calculate 
the lifetimes of d bands in copper fully solving the QP 
Eq. (|^) in the complex plane without relying on any ana- 
lytic continuation. The convergence of the results is care- 
fully checked. Our results are significantly different from 
those obtained within a DFT self-energy based approach 
and are in good agreement with experiments. The care- 
ful analysis of the physics underlying the GW decay of 
quasiholes will shed light into the quantum-mechanical 
mechanism behind the electron dynamics in noble met- 
als. 

In our approach, we start by solving Eq. (|^) on the 
real axis, as it is commonly done in QP bandstructure 
calculations In this way we obtain e^^'", a first 

guess for the real QP energies: 
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where E„k (w) = (nk|i; (r r', w) |nk) and V^'^ = 
(nk|T4;c (r) |nk). Even if Eq. (|l|) in principle requires 
a diagonalization with respect to the band index n, it 
can be reduced to the form of Eq. ^ because the com- 
puted off-diagonal matrix elements of S are at least 
two orders of magnitude smaller than the diagonal ones 
<< |S„,„| Vn,n'). The difference between the 
requested exact quasiparticle energy e^jf and the first 
guess defined in Eq. (||) is due to the fact that the self- 
energy has an imaginary part, namely: 
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Since the quasiparticle concept holds when 
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^nk'^ is a natural start- 
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ing point to get the QP excitation. Note that being Z„k 
complex Eq. will also slightly modify the real part 



of e 
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Thus the corresponding linewidth (r„k) and 



lifetime (r„k) are given by 
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A remarkable property of Eq. (||) is that, being e'^^'^ solu- 
tion of Eq. (||) , only the Re [Znk] is needed to define the 
QP lifetime in contrast with the common expansion of 
S„k {^) around the DFT energy e^w^- is the usual 



renormalization factor, referred to the initial QP guess 
instead of to the DFT eigenvalue. A first approximate 
solution of Eqs. (Q-|2|) can be obtained by fully neglecting 



the QP correction Re 
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and hence Z„k = 1- 



, i.e., by assuming 
The corresponding 



lifetime is T„k ~ {2Im [E„k (e^if "^)] } and is usually 
referred in the literature as the "on mass-shell" GqWo 
lifetime |l^, because the input energies of the QP 
equation are supposed to remain constant, and are subse- 
quently used to calculate the lifetimes. Another approach 
uses Eq. with real Z„k, and e^^'^ corresponding to a 
LMTO-GVF bandstructure This method has con- 
firmed the overestimation of the top d-bands lifetimes 
found in the "on mass-shell" approximation JTsl . How- 
ever, we stress that the approximated real QP energies 
^nif used in Ref. ||^ are not solution of Eq. (Q) . As a 
consequence, Pnk, obtained from Eq. (^), contains and 
additional term proportional to Im [Z„k] that usually is 
neglected ||^ or, it is very small. However for the case 
of noble metals it can be as large as 20 meV for the top 
d-bands (see Fig. |^), accounting for 10-40% of the to- 
tal electronic linewidth. This term reduces the lifetimes, 
in agreement with the experiment and with the results 
presented below. 

Al the present calculations of Green's function and 
screened interaction have been performed using a plane 
wave basis. We have used norm-conserving pseudopo- 
tentials ||l|, [l^ in the DFT-LDA calculation, with 
a 60 Ry energy cutoff, corresponding to ~ 800 plane 
waves |l^. A particular care has been devoted to the 
check of the effect of the Lorentzian broadening r] in- 
cluded in the screening function W [uj) for numerical rea- 
sons (see Ref. |l^ for details). Calculations with differ- 
ent broadening values have shown an increase up to 50% 
of the lifetimes when rj is reduced from 0.1 to O.OOSeV. 
This dramatic numerical effect can be avoided by using 
ri = O.OOSeV, yielding results which coincides with those 
extrapolated at = eV. This is an important point 
since most calculations presented so far have been done 
using dampings of 0.1 eV or larger. 

In Fig. 1^ we present the "on mass-shell" GoWo result 
compared with the full solution of Eqs. (||-^ (simply re- 
ferred as to GqWo), and with the experimental photoe- 
mission results Q. The "on mass-shell" GqWo yields 
lifetimes which are four times too large at the d-bands 
top, and three times too small at the d-bands bottom 
(not shown in Fig. [|). The GqWo results, instead, are 
in rather good agreement with experiments. Actually, 
they are systematically above the experimental values, as 
expected for the contribution of higher-order electron- 
electron and residual phononic and impurity contribu- 
tions. Moreover, as shown in Tab. |, the energy positions 
of the quasiparticle peaks are well reproduced in GoWqj 
while in the "on mass-shell" calculation the eigenvalues 
(coincident with DFT ones) span a larger range of en- 
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ergies, reflecting the known DFT overestimation of the 
hnewidth of rf-bands The origin of the big diff'er- 

ence between the two lifetime calculations stems from the 
large QP self-energy corrections to the d-bands of cop- 
per, completely neglected in "on mass-shell" calculations. 



Wo (r,r'; w') is the dynamically screened potential (con- 
volution of the inverse dielectric function {r,r';uj') 
with the bare Coulomb potential). We can write explic- 
itly Eq. (||) in terms of a summation over the DFT states 
embodied in Go as: 



100 



80 



60 



40 



20 




3.5 
4.0 



> 
E 



E-Ep [eV] 



FIG. 1: Lifetimes of selected d-bands of copper. Diamonds with 
error bars: experimental data from Ref. H. The theoretical results 
obtained in this work are reported for different level of iterations 
within the the GiWo quasiparticle approximation. Full line; GqWo- 
Dotted line: "on mass-shell" GqWo- Dashed line: "on mass-shell" 
GlWo- AH theoretical quasiparticle energy position are aligned to 
the corresponding experimental values. For a comparison of the 
band positions see Tab. hi. The lifetimes are systematically above 
the experimental values, as expected for the contribution of higher- 
order electron-electron and residual phononic contributions. 
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TABLE I: Theoretical bandwidths (in eV) and band energies for 
copper, at high-symmetry points and for various iterations of the 
GiWo quasiparticle approximation. There is a striking agreement 
with the experiment at the Go Wo level [ho||, but this worsens when 
the number of iterations is increased, showing the potential impor- 
tance of including also, vertex corrections. The experimental values 
are taken from ref. |19| . 



The proper inclusion of these non trivial self-energy 
corrections change the electronic decay channels of the d 
levels. In particular, we expect a strong effect on the top- 
most bands. Let us develop this idea further, by looking 
at the GqWo expression for the self-energy: 
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Go(rr',c^')W'o(rr',c^-c^')- (6) 



W^nk^n'(k-q) (^^) = J drdv' ^p^^ (r) ■0„'(k-q) (r) 

VF(rr',c.)V':,(k_q)(r')V'nk(r'), (8) 

being /„k the Fermi occupations and 9 (w) the step func- 
tion. In Eq. (0) a quasiholc in the state |nk) looses energy 
via transitions to all the possible occupied states with 
(higher) energy f^^i^^^^y the energy difference is dissi- 
pated by the the screening cloud (described by W (w)) 
that surrounds the DFT hole |nk}. Thus the GqWq tran- 
sitions contributing to the hole linewidth look like sin- 
gle particle transitions from a quasihole to a DFT hole. 
In higher order of Hedin's equations, i.e. including ver- 
tex corrections, this interpretation of transitions looses 
meaning. This is due to self-energy effects on the Green's 
function of Eq. (|^) and interactions of the screening cloud 
with the state |n' (k — q)} (vertex corrections |l2|) that 
in Eq. (0) are neglected. 

For the top of the occupied d bands the GoW^ cal- 
culation yields negative QP corrections [|o|; this means 
that Eq. (^ contains also contributions coming from the 
decay of the QP d band to exactly the same DFT band. 
We will refer to these transitions as "intraband" decay 
channels. These contributions are important because the 
(i-bands of copper are flat, hence the corresponding den- 
sity of states is large. Moreover the screened interaction 
between d-bands is strong, as the d-states are spatially 
localized and screening is less effective at small distances. 
In the "on mass-shell" calculation the states |n' (k — q)) 
appearing in Eq. (^ and the quasiparticle states corre- 
spond to the same DFT eigenvalues; this means that the 
"intraband" decay channels occur at zero energy, where 
the low-energy Drude tail of the dielectric function dom- 
inates (and hence the screened interaction is vanishing). 
This leads to the usual interpretation of the long calcu- 
lated lifetimes at the d bands top as due to the fact that 
these d states can only decay to s/p states. These matrix 
elements are smaller than those involving states with the 
same Z-character. 

Even if the results of the Go Wo calculation are in 
rather good agreement with the experiments, a natural 
question about the physical meaning of the "intraband" 
decay channels arises. Being not at self-consistency, the 
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system is described within GqWq on the basis of quasi- 
particle states and DFT states (those involved in the hole 
decay) with different energies. To remove this inconsis- 
tency Eq. (|^) should be solved iteratively, until converged 
QP energies are obtained. However, as we show below, 
iterations beyond the usual GqWo level worsen both the 
imaginary and real parts of the calculated QP energies. 

So far fully self-consistent GW calculations have been 
performed only for the homogeneous electron gas |Q and 
for simple semiconductors and metals pH , yielding worse 
spectral properties than those obtained in the non self- 
consistent GqWq. The construction of a self-consistent 
GW self-energy is a formidable task even for the simple 
systems mentioned above. In copper, already the update 
of the screening function is rather demanding, due to the 
presence of localized d orbitals that imply a large cut-off 
in the plane wave expansion. To test the effect of self- 
consistency on the QP energies and lifetimes we use a 
simplified GiWo method, where the self-energy operator 
is defined as 

S.(rr',c.)^ j'-^Gdvv',uj')W(o)ijv\oj~oj'), (9) 

being i is the iteration number. Gt involves the QP ener- 
gies obtained from without considering renormal- 
ization factors, lifetimes, and energy structures beyond 
QP peaks. As the QP bandstructure resulting from the 
first iteration is already in excellent agreement with ex- 
periment, our next step is to perform an "on mass-shell" 
GiWq calculation. The resulting lifetimes are compared 
with experiment in Fig.|l|. One sees that forcing QP en- 
ergies to appear also in the states involved in Eq. i.e. 
the states forward which the quasihole decays, yields life- 
time results similar to those of the "on mass-shell" GqWq 
method. The same overestimation of the lifetimes of the 
top of the d bands is observed, confirming that a good 
agreement with experiment depends on the inclusion of 
the "intraband" decay channels described by Eq. (0) . 

A further question could now be addressed and is how 
self-consistency affects the quasiparticle bandstructure 
obtained within GoWq. As shown in Tab. |, at higher 
iteration orders of the quasiparticle GW equation the re- 
sulting energies worsen. The d-bands width decreases, 
reducing the agreement with experiment. 

These results show that GqWq describes correctly the 
rf-holes lifetimes as far as "intraband" decay channels be- 
tween (i-like states are included. Those transitions, how- 
ever, imply an inconsistency between the quasiparticle 
initial states and the DFT final states of the hole decay, 
as shown in Eq. (0). A self-consistent solution of Dyson 
equation removes this inconsistency, worsening, however, 
the agreement with the experimental results. This is a 
clear indication of the need of including vertex correc- 
tions in the self-consistent procedure. As shown for sim- 
ple metals ||2^ , vertex corrections would partially cancel 
the dressing of the Gq Green's function of Eq. (^), restor- 



ing the "intraband" decay channels and, consequently, 
the Go Wo results. 

In conclusion, we have shown that the lifetimes of d- 
holes in copper, calculated within the Go Wo method, are 
in good agreement with the experimental results, and can 
be obtained within the very same scheme which yields a 
good quasiparticle bandstructure. In contrast, further it- 
erations of the QP equation beyond the Go Wo level yield 
worse results, for both the real and imaginary parts of 
self-energy. This can be explained by the need of includ- 
ing also vertex corrections together with self-consistency. 
This result is quite general and should apply to all met- 
als in which one has two or more sets of electronic states 
with different degrees of spatial localization. 
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